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We study the time evolution of velocity and pressure gradients in isotropic turbulence, by quan- 
tifying their decorrelation time scales as one follows fluid particles in the flow. The Lagrangian 
analysis uses data in a public database generated using direct numerical simulation of the Naiver- 
Stokes equations, at a Reynolds number Re\ « 430. It is confirmed that when averaging over 
the entire domain, correlation functions decay on timescales on the order of the mean Kolmogorov 
turnover time scale, computed from the globally averaged rate of dissipation and viscosity. How- 
ever, when performing the analysis in different subregions of the flow, turbulence intermittency 
leads to large spatial variability in the decay time scales. Remarkably, excellent collapse of the auto- 
correlation functions is recovered when using the 'local Kolmogorov time-scale' defined using the 
locally averaged, rather than the global, dissipation-rate. This provides new evidence for the validity 
of Kolmogorov's Refined Similarity Hypothesis, but from a Lagrangian viewpoint that provides a 
natural frame to describe the dynamical time evolution of turbulence. 

PACS numbers; 02.60.Cb,02.70.Hm,47.27.Ak,47.27.ek,47.27.Gs 



Understanding the universal features of the dynamics 
of turbulence [l|, 3 continues to be a formidable prob- 
lem in classical and statistical physics. The original 
Kolmogorov 1941 theory relied on the overall averaged 
dissipation-rate (e) to predict, among others, the scaling 
properties of the energy spectrum. It was extended to 
account for intermittency by the introduction of the re- 
fined Kolmogorov similarity hypothesis (RKSH) In 
this extension, attention is placed on conditional statis- 
tics based on the locally averaged dissipation-rate in some 
particular subregion of the flow, such as a sphere or box 
of size r. This local dissipation-rate, usually denoted by 
Cr, is defined according to 
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where V is the volume of the subregion TZr{x) of size 
r centered at x, ly is the kinematic viscosity of the 
fluid, and Sij is the strain-rate tensor defined as Sij — 
{dui/dxj + duj/dxi)/2 (where Ui is the velocity field). 
The longitudinal velocity increment at scale r is defined 
as 5rU = [ui(x + r) — Uj(x)](ri/r), and the RKSH states 
that in the inertial range of turbulence the statistics of 
5rU depend on r and so that from dimensional anal- 
ysis moments of 5rU, conditioned upon a fixed value of 
Cr will scale as {5rU^\er) = Cp{rerY^^^ according to Kol- 
mogorov's 1941 postulate. Anomalous scaling then re- 
sults from the additional global averaging and anomalous 
scaling behavior of moments of e^- 

Extensive literature to test RKSH has focused mainly 
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on velocity increments 0, 0, H, 0, H, 01 or acceleration 
p^ . In recent years there has been growing attention 
placed in the dynamical evolution of the velocity gradi- 
ent tensor A {Aij = dui/dxj) due to the fact that A 
provides rich information about the topological and sta- 
tistical properties of small-scale structure in turbulence. 
The Lagrangian time evolution of A can be obtained by 
taking gradient of the NS equation pd| : 
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where d/dt stands for Lagrangian material derivative and 
p is the pressure divided by the density of the fluid. The 
flrst term on the right-hand side of Eq. ([2]) denotes the 
nonlinear self-interaction of A, the second term is a ten- 
sor called pressure Hessian Pij = d^p/dxidxj, and the 
third is the viscous term. Assuming the pressure Hessian 
is isotropic (i.e. neglecting df^p — dlj^p Sij/3) and neglect- 
ing the viscous term lead to a closed formulation for A, 
the so-called Restricted-Euler (RE) equation, which has 
analytical solutions for the full tensor-level time history. 
Remarkably this simple system is already sufficient to 
explain a number of non-trivial geometrical trends found 
in real turbulence [H!, [l^l- Nevertheless, the RE system 
leads to nonphysical finite-time singularities because the 
self-stretching is not constrained by any energy exchange 
or loss mechanism in the system. To develop models 
for such energy exchange mechanisms, there is a need to 
better understand the time evolution of A, as one follows 
fluid particles across a turbulent flow. As will be shown 
in this paper, the RKSH can play a crucial role in deter- 
mining the characteristic time scales of this evolution in 
different parts of the flow. 

Much effort at regularizing the RE system to avoid the 
nonphysical singularity has been made [l^, 0, [l^, [l^, [13, 
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llSj l in the past two decades. Despite progress, the inabil- 
ity to fully account for the anisotropic pressure Hessian 
and viscous effects continues to limit the accuracy of the 
exis ting models for the evolution of velocity-gradient ten- 
sor Il9|. In various models of the Lagrangian dynamics 
the characteristic correlation times along 
Lagrangian trajectories and their scaling with Reynolds 
number play a central role. For instance in the model 
based on the 'recent fluid deformation closure' the 
Lagrangian pressure Hessian tensor is assumed isotropic, 
based on the idea that any causal relationship between 
initial and present orientations will be lost after a char- 
acteristic Lagrangian correlation time-scale of the ten- 
sor A. The usual expectation is that the characteristic 
correlation time-scale of A is the Kolmogorov time-scale 

Since the dynamical equation for the veloci ty g radient 
tensor A (Eq. ^ and the existing models [l^, [UlU, [H, 
[TtI . [Tsj are written for the full tensor, it is of interest to 
quantify the temporal correlation function of each tensor 
element but to do so in a fashion that is coordinate sys- 
tem invariant. We thus use the tensor-based Lagrangian 
time correlation function [23j of a second-rank tensor C, 
defined as 
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where r is the time-lag along Lagrangian trajectories and 
(• • •) may represent ensemble or volume averaging for 
homogeneous turbulence. Here tensor elements arc as- 
sumed to have zero mean for isotropic turbulence. For 
a statistically steady-state process, when pc{t) does not 
depend upon t^, averaging can also be done over tg. The 
correlation function pc (''") provides a frame- invariant de- 
scription of the auto-correlation structure of tensor ele- 
ments, appropriately summed over all directions. 

The fluid acceleration is another variable of great inter- 
est. The Lagrangian velocity and acceleration statistics 
have been investigated both experimentally [2^ and nu- 
merically (25I |26| . The acceleration is associated more 
closely with inertial-rangc and small-scale structures. It 
has been established that the acceleration is dominated 
by the pressure gradient Vp and viscous forces are negli- 
gible away from boundaries [23,[11|- Therefore, we study 
the time correlations associated with the pressure gradi- 
ent. Specifically, similar to the tensor correlation func- 
tion defined in Eq. ([3]), for the pressure gradient correla- 
tion pvp{t) we use the inner product of the vector at to 
and ^0 + 

We measure Pc{t) for C = A, S or 



Aji)/2) and pvpi^) using data 



from pseudo-spectral direct numerical simulation of 
forced isotropic turbulence on a 1024'^-node peri- 
odic domain, with Re\ = 433. The 27 Terabytes 
of data are stored in a public database format (see 
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the trajectories, such as velocity and pressure gradients 
through the second-order Runge-Kutta particle-tracking 
algorithm [2^. The interpolation of required velocities 
uses Sth-order Lagrange polynomials in space and 
piecewise cubic Hermite polynomial interpolation in 
time, as implemented in the predefined functions in the 
database j3(j |. 

The Lagrangian time correlations for S and Q are 
shown as open squares in Fig. [ija) (the meaning of 
the other lines are explained below). As it can be 
seen, the strain-rate auto-correlation decays fairly rapidly 
and reaches almost zero near 5-6 Kolmogorov timescales, 
much faster than that of the rotation-rate [231 . 




and can be accessed using web-service tools. We track 
fluid particles and extract Lagrangian information along 



FIG. 1: Lagrangian auto-correlations of strain- and rotation- 
rate tensors computed from Eq. (|3]). Open squares are for 
global average over randomly located particles, whereas dif- 
ferent lines correspond to subregions of the flow characterized 
by certain er. Time-lag is normalized using (a) the global 
Kolmogorov time scale tk and (b) the local time-scale TK,r- 



It is known [3| that turbulence is highly intermittent, 
with regions displaying strong fluctuations, interspersed 
with less turbulent regions. In order to study intermit- 
tency from the viewpoint of the Lagrangian time evo- 
lution, we compute the correlation functions based on 
fluid particles that originate from various subregions of 
the flow. The subregions are characterized by the lo- 
cal dissipation-rate e^. We thus define the conditional 
time-correlation functions based on e^- In Eq.Q. the 
three global averages are replaced by conditional aver- 
ages, e.g. {Cij{to)Cij{to + t)) {Cij{to)Cij{to + T)\er). 
In these averages the initial position of particles con- 
tributing to the average at time Iq are sampled from 
several local boxes of size r that have a prescribed lo- 
cally averaged dissipation-rate Cr (in practice a range of 
values in a prescribed 'bin' is used). We consider four 
scales r in the inertial range of turbulence correspond- 
ing to 16, 32, 64 and 128 grid-points of the DNS or 
r = SArjK, GSrjK, 13677;^ and 276riK respectively where 
rjK is the Kolmogorov length scale. Four bins of val- 
ues are chosen, centered at: Cr =0.0076, 0.0186, 0.0924 
and 0.148 for 16-cubes, 0.0096, 0.0186, 0.092 and 0.148 
for 32-cubes, and 0.0338, 0.0519, 0.0936 and 0.122 for 
128-cubes. For scale r = 136t]k (64-cube), we have five 
bins which are centered at = 0.0371, 0.0458, 0.0659, 
0.0859, and 0.1487. The overall mean dissipation-rate 
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over the entire dataset is (e) = 0.093 jSOj . 

To illustrate the situation, in Fig. [5] we show 12 repre- 
sentative 64-cubes placed inside the 1024'^ domain, with 
50 sample fluid particle trajectories emanating from each 
and progressing during a time equal to 27 tk- The re- 
quired averages are taken over all these trajectories as 
well as over several cubes for which is in a bin's pre- 
scribed range. The PDF of the dissipation-rate is shown 
in the insert in Fig. [2] together with the bins or values 
used for the conditional averaging. 



log,„(e/<E>) 




FIG. 2: Sample particle trajectories starting from 12 ran- 
domly selected 64-cubes characterized by local dissipation- 
rate at the initial time. Contours on the background planes 
is local (normalized) dissipation-rate (at the initial time) in 
logarithmic units showing its intermittent but structured dis- 
tribution at the smallest scales of the flow. The lower right 
insert shows the PDF of the dissipation-rate together with 
the bins or tr values used for the conditional averaging. 



For the case of r = 13677^ (64-cube), the various 
lines in Figure[T]Ja) show the conditional auto-correlation 
functions so computed, plotted as function of time-delay 
scaled by global Kolmogorov time scale tk ■ Averages are 
evaluated over 6 x 10^ particles in each bin. Tests varying 
the number of particles show that the auto-correlation 
functions are well converged with this number of parti- 
cles. There are noticeable differences in the results de- 
pending on e^, for both the S and il auto-correlation 
bmctions. For larger values of e^, i.e. in regions of more 
intense turbulence activity, the decay rate of the auto- 
correlation functions is faster. This shows that the global 
Kolmogorov time-scale tk = {y/{e)Y/'^ does not deter- 
mine the local evolution of patches of turbulence, even 
in a statistical sense when conditional averaging is used 
that segregates different types of fluid regions. 

As a next step, consistent with the RKSH, we de- 
fine a 'local Kolmogorov time-scale' based on the local 
dissipation-rate according to 

TK,r = {v/erY/'^. 



A Lagrangian version of the RKSH would state that 
the temporal auto-correlation functions should be a uni- 
versal function of a time-delay r normalized by a 'local 
Kolmogorov time scale' TK,r- Results shown in Fig. [T^b) 
thus scaled by the local time-scale for each tensor show 
excellent collapse. In this plot we also show the auto- 
correlation for the velocity-gradient tensor A which falls 
between the results of its symmetric and antisymmet- 
ric parts. Again, the auto-correlation of the rotation- 
rate tensor decays much more slowly than that of the 
strain-rate tensor, which can be attributed to the ef- 
fect of the deviatoric part of the pressure Hessian as 
discussed more in depth in Ref. [23|. The deviatoric 
part of the pressure Hessian causes non-local effects on 
the strain-rate evolution but has no effects on the evolu- 
tion of rotation-rate. Figure [3] shows similar results for 
the pressure gradient. Also here, the collapse is much 
more improved when using the local time-scale that cor- 
responds to the local dissipation-rate, showing that the 
Lagrangian KRSH works for both the velocity and pres- 
sure gradients, the latter corresponding closely to the 
fluid acceleration [13, [H]. 

We point out that in recent work [s^ the Lagrangian 
KRSH has also been shown to hold in the context of mo- 
ments of two-time velocity increments. In their analysis, 
the authors use a rate of dissipation e,- averaged over 
temporal domains of duration r along the particle trajec- 
tory. The 'fully Lagrangian' quantity e,- averages, and 
thus connects, dissipation at various times. From a dif- 
ferent, perhaps complementary point of view, the present 
analysis based on the more traditional spatial average of 
dissipation at a single (initial-condition) time facilitates 
interpreting the results in terms of 'causality' based on a 
fixed initial condition. 
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FIG. 3: Lagrangian auto-correlations of pressure gradient, 
time-lag is normalized using (a) global time-scale tk and (b) 
local time scale tk r- 



Next, we explore the validity of the Lagrangian KRSH 
for different box sizes r. Figure 2] shows the auto- 
correlations of strain-rate tensor ps for the four box sizes 
r and described above. Now the scatter is even more 
pronounced in (a) where time-lag is normalized by the 
global averaged Kolmogorov time scale tk, mainly due 
to the fact that the range of values of for smaller values 
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of r is larger due to intcrmittency. Once again, excellent 
collapse can be observed in (b) where the time-lag is nor- 
malized by the local Kolmogorov time scales TK,r, even 
though the curves correspond to different r and e,. values. 
This is confirmed further by comparing the results with 
similar local Cr values but different r: e.g. the first two 
curves from the right to left shown in Fig. 3] (a) corre- 
spond to = 0.0076 for a f 6-cubc and = 0.0096 fora 
32-cube and they collapse quite well. Similar results (not 
shown) are obtained for rotation rate tensor. 




Present results provide a new and more dynamical in- 
terpretation of the Kolmogorov Similarity Hypothesis: 
when focussing on particular subregions of turbulence, 
the dynamics proceed according to time-scales dictated 
by the locally averaged rate of dissipation. The classical 
view of the KRSH focussed on the moments of velocity 
increments in turbulent fields that had already developed 
due to past dynamics. Here we have shown that the dy- 
namical time evolution of turbulence, its rate of change 
and associated decorrelation time scale, depends on the 
local rate of dissipation, averaged over a volume com- 
parable to the volume containing the sample of initial 
particle locations. Results also suggest that Lagrangian 
models (e.g. for the velocity gradient tensor 17| or the 
acceleration 
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that require specification of a charac- 
teristic time scale should not use the global, but the local 
time scale. 



FIG. 4: Lagrangian auto-correlations of strain-rate tensor 
computed from Eq. Q for four box-size values in the inertia! 
range r, and for four values (except for 64-cube which has 5 
bins of Er). The lines are solid lines : 16-cube; • • • : 32-cube; 

: 64-cube; — • ■ : 128-cube; open squares: global average. 

In (a) time is scaled using global time-scale, tk, in (b) using 
the local time scale, tk.t- 
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